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Abstract 



We study the soft limit of one-loop QCD amplitudes and we derive the process- 
independent factorization formula that controls the singular behaviour in this 
limit. This is obtained from the customary eikonal factorization formula valid 
at tree (classical) level by introducing a generalized soft-gluon current that 
embodies the quantum corrections. We compute the explicit expression of the 
soft-gluon current at one-loop order. It contains purely non-abelian correlations 
between the colour charges of each pair of hard-momentum partons in the 
matrix element. This leads to colour correlations between (two and) three 
hard partons in the matrix element squared. Exploiting colour conservation, 
we recover QED-like factorization for the square of the matrix elements with 
two and three hard partons. 
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1 Introduction 



Higher-order computations in perturbative QCD can be performed by using three main 
tools: exact calculations at a fixed order in the QCD coupling as = S's/^^r, analytic re- 
summed calculations and parton shower event generators. At present, the accuracy of these 
tools is respectively limited to the next-to-leading order (NLO), to next-to-leading loga- 
rithmic (NLL) accuracy and to the dominant soft and collinear enhanced contributions. 
An extensive and updated list of references can be found in Refs. [jl], 0] . 

Apart from important conceptual and technical differences, these tools are based on 
a common ingredient: the universal factorization properties of QCD amplitudes in the 
infrared (soft and collinear) region. The lowest-order version of the soft and collinear 
factorization formulae P, ^ was indeed exploited to develop these tools to their present 
theoretical accuracy. Higher-order versions of the factorization formulae are required to 
progress towards next-to-next-to-leading order (NNLO) calculations, resummation of next- 
to-next-to-leading logarithmic (NNLL) terms and inclusion of subdominant contributions 
in parton showers. 

In recent years several groups have contributed to extending infrared factorization to 
higher perturbative orders fl^T6|. 

The general factorization properties of tree-level and loop amplitudes in the limit where 
two or more partons become collinear were studied in Refs. [|]-|^. At tree level, the singular 
factors for the collinear splitting of one parton into three were explicitly computed in 
Refs. 1^, |10|- The splitting function corresponding to the clustering of four collinear 
gluons is also known [|TT|]. The one- loop kernels for the collinear splitting of one parton into 
two were obtained in Refs. |1I2|-|T4|. 



The mixed soft-collinear limit |^ can be studied by exploiting the coherence prop- 
erties of QCD radiation [^]. Using QCD coherence, the singular behaviour in the soft- 
collinear limit can be treated [j^ by combining the singular factors that separately control 
the collinear and soft limits. 

The soft limit is physically more involved than the collinear limit. Long- wavelength 
(soft) gluons can spread the colour fiow over large distances, thus leading to (non-local) 
colour correlations. Kinematics and colour factors turn out to be deeply entangled in the 
soft-factorization formulae. 

The tree-level factorization formulae for the emission of two soft gluons were indepen- 
dently derived in Refs. |15] and 0. The soft (and collinear) singular behaviour of two-loop 



amplitudes was studied in Ref. ||T6[. 

In this paper we consider the limit in which a soft gluon is radiated from one-loop 
amplitudes. The limit was first investigated in Refs. [|, 14 1. The formalism used by 



this group is based on the decomposition of the one-loop matrix elements in colour sub- 
amplitudes |T^. The colour-subamplitude decomposition depends on the type of external 



partons, and these authors derived the explicit expressions of the one-loop soft-gluon con- 
tribution to colour subamplitudes with m external gluons |T2| and with m external gluons 



plus a qq pair [14 . 
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In the present paper the soft hmit is studied by means of a completely independent 
and general method. We apply the eikonal approximation and soft-gluon insertion rules 
to perform infrared factorization directly in colour space. Within this formalism, the soft 
limit of tree-level amplitudes is described by a factorization formula written in terms of 
a soft-gluon (or eikonal) current (see Ref. @] and Sect. ||) that describes colour radiation 
in the classical approximation. We show that the factorization formula can be extended 
to loop amplitudes by introducing a generalized soft-gluon current that embodies (non- 
abelian) quantum corrections. As at tree level, the soft-gluon current only depends on the 
colour charges and momenta of the external partons in the loop amplitude. We compute 
the explicit expression of the soft-gluon current at one-loop order. 



Our results for the soft limit of one-loop amplitudes agree with those in Refs. |T^, |14 
for the particular cases considered therein. The one-loop gluon current that we obtain is 
derived and presented in a general and process-independent way; it can be applied to any 
one-loop amplitude. In particular, we can easily show that colour and kinematic factors can 
be completely disentangled in the computation of the soft limit of the square of one-loop 
matrix elements with two or three external QCD partons. This simplified factorization 
structure is particularly useful for the NNLO calculation of 2-jet and 3-jet cross sections in 
e~^e~ annihilation. 

Infrared-factorization properties at the lowest perturbative order have been known for 
a long time [Q. Combining the results obtained here with those in Refs. []6|-[16|, the general 
structure of the infrared singularities at the next perturbative order is also completely and 
explicitly known. It can be used to improve the accuracy of perturbative QCD calculations. 

The paper is organized as follows. In Sect. |^ we recall the known results for the soft 
limit of tree-level amplitudes. In Sect. || we present the process-independent factorization 
formula at higher perturbative orders and we discuss in detail its features at one-loop order. 
In Sect. § we prove the factorization formula at one- loop order and we derive the explicit 
expression of the one- loop soft current. In Sect. ^ we apply our results to the squared 
amplitudes of processes with two or three hard partons. 



2 Soft-gluon factorization at tree level 

We consider a generic scattering process that involves m external QCD partons {massless 
quarks and gluons) with momenta pi, . . . , pm and an arbitrary number and type of particles 
with no colour (photons, leptons, vector bosons, ...). Note that, by definition, we always 
consider incoming and outgoing parton momenta in the physical region, i.e. any pi is 
massless, with positive-definite energy (in particular, Pi-pj > 0). The corresponding matrix 
element is denoted by M.{pi, . . . ,Pm) and the dependence on the momenta and quantum 
numbers of non-QCD particles is always understood. 

The matrix element has the following loop expansion: 

M{pu...,Pm) =M^^\pi,...,Pm)+M^^\pi,...,Pm)+... , (1) 

where denotes the tree-level contribution, A4'^^^ denotes the one-loop contribution, and 
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the dots stand for higher-loop corrections. Note that we always consider unrenormalized 
matrix elements. Thus Eq. (|I]) is a power series expansion in the bare QCD coupling gs 
and, in particular, Ai^^^ is the unrenormalized one-loop amplitude. 

We simultaneously regularize ultraviolet and infrared singularities by using dimensional 
regularization. Apart from introducing the dimensional-regularization scale n through the 
replacement gs — > (^s A*^, the key ingredient of dimensional regularization is the analytic 
continuation of loop momenta to d = A — 2e space-time dimensions. Having done this, we 
are left with some freedom regarding the dimensionality of the momenta of the external 
particles as well as the number of polarizations of both external and internal particles. 
This leads to different regularization schemes within the dimensional-regularization 

prescription. The regularization- scheme dependence of one-loop amplitudes was studied in 
detail in Refs. [^, |2^. All the results on the soft behaviour presented in this paper do not 
explicitly depend on the dimensional-regularization scheme (see Sect, ^for a brief discussion 
of different regularization schemes): the scheme dependence is implicitly embodied in the 
expressions of the tree- level and one-loop matrix elements Ai^'^^ and A4^^\ 

The emission of a soft gluon does not affect the momenta and spins of the radiating 
hard partons. However, it does affect their colour because the gluon always carries away 
some colour charge, no matter how soft it is. Unlike the case of soft-photon emission in 
QED, soft-gluon emission thus does not factorize exactly and leads to colour correlations. 

To take into account the colour structure without referring to any particular choice of 
basis colour vectors (such as, for instance, the decomposition in colour subamplitudes [ll4|), 
we use a general notation (see e.g. Ref. [^). The dependence of the matrix element on 



the colour indices ci, . . . , of the QCD partons is written as 

Mc,,...,cAPi,. . . ,Pm) = {Ci,...,Cm\M{pi,...,Pm)) ■ (2) 

Thus { |ci, . . . , Cm) } is an abstract basis in colour space and the ket \J^{pi, . . . ,pm)) is a 
vector in this space. According to this notation, the matrix element squared |A^p (summed 
over the colours and spins of the partons) can be written as 

\M{pi,...,Pm)\^ = {M{pi,...,Pm)\M{pi,...,Pm)) ■ (3) 



To describe the colour correlations produced by soft-gluon emission, we associate a 
colour charge T, with the emission of a gluon from each parton i. If the emitted gluon has 
colour index a [a = 1, ~" 1); the colour-charge operator is: 

T, ^ {a\ Tt (4) 

and its action onto the colour space is defined by 

(a, Ci, . . . , Cj, . . . , Cm\Ti\bi, . . . , 6j, . . . , bm) = Scibi---T^ibi---^Cmbm 5 (5) 

where T^!^ = ifcat (colour-charge matrix in the adjoint representation) if the emitting parton 
z is a gluon and T^^ = (colour-charge matrix in the fundamental representation with 
a, P = 1, . . . , Nc) if the emitting particle z is a final-state quark or an initial-state antiquark 
(T^fl = = — tJ^Q,, in the case of a final-state antiquark or an initial-state quark). 
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The colour-charge algebra is|^ 

T-T^ = T,.T,=T,-T, ifz^j; = (6) 

where Cj is the Casimir operator, i.e. Ci = Ca = ii i is a. gluon and Cj = Cp = 
{N^ — l)/2Nc if i is a quark or antiquark. 

Note that, by definition, each vector \J^{pi, . . . ,Pm)) is a colour-singlet state. Therefore 
colour conservation is simply 

m 

J2T^\M{pu...,pJ) = 0. (7) 

i=l 

We can now recall the behaviour of the tree-level matrix element Ai^^\q,pi, . . . ,Pm) in 
the limit where the momentum q of the gluon becomes soft. Denoting by a and e^{q) the 
colour and the polarization vector of the soft gluon, the matrix element fulfils the following 
factorization formula |^ 

{a I M^'\q,pr, . . . ^ gsf^'e'iq) Jl^'\q) \ M^'\p,, . . . , (8) 

where \ AA^^\pi^ . . . ,Pm)) is obtained from the original matrix element by simply removing 
the soft gluon q. The factor jf^{q) is the tree-level soft-gluon current 

which depends on the momenta and colour charges of the hard partons in the matrix 
element on the right-hand side of Eq. (^. The symbol ' ~ ' means that on the right-hand 
side we have neglected contributions that are less singular than 1/g in the soft limit g ^ 0. 
Note that Eq. is valid in any number d = 4 — 2e of space-time dimensions, and the sole 
dependence on d is in the overall factor n^. 

The factorization formula @ can be derived in a simple way by working in a physical 
gauge and using the following soft-gluon insertion rules. The coupling of the gluon to any 
internal (i.e. highly off-shell) parton in the amplitude Ai^^\q,pi, . . . ,Pm) is not singular 
in the soft limit; it can thus be neglected. The soft-gluon coupling to any external or, in 
general, nearly on-shell parton with colour charge T and momentum p can be factorized by 
using the eikonal approximation, that is by extracting the contribution gsfJ^'^^p^T for the 
vertex and the contribution l/{p + q)"^ — + 2p ■ q) for the propagator. Note that the 
eikonal vertex only depends on the momentum and colour charge of the radiating parton: 
it does not depend on either the soft momentum or the spin of the parton. Using the 
eikonal propagator simply amounts to neglecting the terms that are quadratic in the soft 
momentum. 

An important property of the soft-gluon current is current conservation. Multiplying 
Eq. by g^, we obtain 

m 

q''J^;\q) = Y.T^, (10) 

i=l 

§More details on the colour algebra and useful colour-matrix relations can be found in Appendix A of 
Ref. 11. 
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and thus 



g'^jW(g)| ...,pJ) = J2T^\M^'^\p,,...,Pm)) = , (11) 



■t=i 



where the last equahty follows from colour conservation as in Eq. (0). Although the factor- 
ization formula (§) is most easily derived by working in a physical gauge, the conservation 
of the soft-gluon current implies that Eq. (H) is actually gauge invariant. Any gauge trans- 
formation is equivalent to an addition of a longitudinal component to the polarization 
vector of the soft gluon through the replacement e'^{q) —>■ e^{q) + \q^. Nonetheless the 
factorization formula (|^) is invariant under this replacement, because of Eq. (|n|). 

Squaring Eq. (H) and summing over the gluon polarizations leads to the well-known 
soft-gluon factorization formula at 0{g1) for the squared amplitude |Q: 

m 

|A^W(g,Pi,...,p„)P^-^7i/.2^2 S.M) |-MSS)(pi,...,P^)r , (12) 

where the eikonal function Sij{q) can be written in terms of two-particle sub-energies Sij = 
{Pi + Pi)^ &s follows 

S^M = wr^^ V = ■ (13) 



2(Pi ■ q) {Pj ■ q) 



Sin S 



iq Jjq 



The colour correlations produced at tree level by the emission of a soft gluon are taken into 

,(0) I 



account by the square of the colour-correlated tree-amplitude on the right-hand 



side. This is given by 

|A^SS)(pi, . . . ,Pm)? = {M^'\p,, . . .,pm) I T, ■ T, I M^'\pu . . . ,p„,)) (14) 



(0) 



3 Soft-gluon factorization at one loop 

The factorized structure of the soft limit of the QCD amplitudes at tree level can be 
generalized to higher loops. The soft limit of the one-loop amplitudes is studied in detail 
in Sect. In this section we anticipate and discuss the final results. 

Our analysis is consistent with the following factorization formula 

{a\Miq,p^,...,Pm))^s^{q)r^{q,e)\M{p,,...,Pm)) [l + 0{gl)] , (15) 

where the symbol ' ~ ' has the same meaning as in Eq. (P) . The matrix element on the 
right-hand side is the all-loop amplitude in Eq. (P and the singular dependence on q is 
embodied in the (unrenormalized) soft-gluon current J^{q,e), which can be expanded in 
loop contributions, i.e. in powers of g^: 

j;(g, e) = /i^ [j; (q) + qI ^^''J^ '•'Hq, e) + . . .] . (16) 
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The term J^^^\q) is the tree-level current in Eq. (P), the term J^^^\q,e) is its one-loop 
correction, and so forth. 

The soft current contains the entire singular dependence in the soft limit and no further 
approximation is performed in Eq. ([T5|) . By this we mean that J^{q,e) behaves as 1/g 
(such as at tree level) modulo any possible enhancement proportional to powers of Inq 
coming from the loop contributions. In particular, the e-dependence of the current can be 
evaluated exactly without performing any e-expansion and, thus, by keeping all the powers 
of Ing coming from higher-loop contributions of the type [qY = 1 + elng + ■ ■ ■ . 

The discussion in Sect. ^ suggests that the factorization formula (^) is valid to any loop 
order. Nonetheless, we have included the term 0{g^) on the right-hand side of Eq. (|T5|) to 
indicate that our explicit proof and calculation do not extend beyond the one-loop order. 

we can obtain 



Expanding both sides of Eq. ([15| ) to one-loop accuracy and using Eq. 
the factorization formula for the soft limit of the one-loop amplitudes: 



M^^\q,Pi,...,Pm)) 



9sf^'e'^iq)[r^^'\q)\M^'\p,,...,Pn)) 

+ 9lf^''j;^'\q,e)\M^'\pu. 



,Pn))] . (17) 



The explicit expression for the one-loop current is 
1 1 r3(l-e)r2(l + e) 



167r2 e2 



r(l - 2e) 

Vi 



Pi-q pj-q 



An Pi -pjC' 



2(Pj ■ q) {Vj ■ <?) 6 ^^''''^ e 



(18) 



We remind the reader that in our notation all the incoming and outgoing momenta are in 
the physical region (any pi has positive-definite energy and pi ■ pj > 0). Thus the complex 
factors e~^'^^^'^ (Aab = +1 if v4 and B are both incoming or outgoing, and Xab = 
otherwise) in Eq. (p!8| ) are the unitarity phases related to the analytic continuation from 
unphysical to physical momenta. 

The result in Eq. (|T8|) explicitly shows that the soft limit of the one-loop amplitudes is 
process- independent, meaning that it does not depend on the momentum and colour flows 
of the internal partons (including the parton circulating in the loop) in the matrix element. 
This simple structure has several interesting features that we comment below. 

The one-loop soft current is proportional to the structure constants fate of the gauge 
group, and thus it is purely non-ahelian. This is in agreement with the absence of higher- 
loop corrections to the soft current in massless QED [ p5| . 

Since Eq. (|TBp is proportional to the factor {Pi /Pi-q—Pj /Pj-q), the one-loop contribution 
to the soft current is conserved: 



. 



(19) 



Combined with the analogous property at tree level, this guarantees that the soft-gluon 
factorization formula ([T5| ) is manifestly gauge-invariant. 
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The double pole in Eq. ( |18|) is the infrared singularity produced by a soft and 
coUinear virtual gluon. To double-pole accuracy, we can use colour conservation (see 
Eq. (0)) to show that the one-loop current is simply proportional to the tree-level cur- 
rent: 



J"«(g,e)|A^(g,Pi,...,p„)) 



167r2 



O 



M{q,pu...,pm)) . (20) 



This behaviour is consistent with the known singularity structure of the one-loop amplitudes 



Beyond the double-pole approximation, the one-loop current contains ||2^ two-particle 
colour correlations. The correlations are induced by the last factor on the right-hand 
side of Eq. (|T^). This factor fully embodies the logarithmic dependence on the soft-gluon 
momentum q and has a simple kinematic interpretation, being related to the transverse 
component q±^ij of the gluon momentum with respect to the longitudinal direction singled 
out by the momenta Pi and Pj of the colour-correlated hard partons: 



2 ^ 2(p^ ■ q) {pj ■ q) 



(21) 



The derivation of the result in Eq. (|1^) (see Sect. H) suggests that multiparticle colour 
correlations will appear in higher- loop contributions to the soft-gluon current (|l^). For 
instance, at two-loop order the soft current contains colour correlations of the type 



£ r rpb rpc rpd 

Jabe Jecd J- i J- j J- k 



(22) 



between three different hard partons i,j and k. 



The soft behaviour of the one-loop amplitudes was first investigated in Refs. ||, |T2|, |T4| 
by using the colour-subamplitude formalism. The explicit expressions of the one-loop soft- 
gluon contribution to colour subamplitudes, with m external gluons and with m external 
gluons plus a qq pair, were derived in Refs. []T2| and |]T^, respectively. The reader can 
straightforwardly check that the result in Eq. ([T8|) agrees with those in Refs. [0, |T^ for 
the particular cases considered therein. The general result in Eq. (|T8|) shows that, although 
colour and kinematics are deeply entangled in the soft region, the soft limit of the one-loop 
amplitudes can be factorized in colour-space in a way that is both (relatively) simple and 
process- independent (in particular, independent of the flavour of the external partons). 
In particular, this general structure is quite useful to show (see below and Sect. H) that 
colour and kinematics can be completely disentangled in the computation of the soft limit 
of the square of one-loop matrix elements with two or three external QCD partons. The 
factorization formula (|T^) can be used to compute the soft limit of the one-loop contribution 
to the square of the matrix element^] Ai{q, {p})- Summing over the polarizations of the 

^In the following the dependence of the matrix element on the momenta pi, . . . , Pm of the hard partons 
is denoted by {p}. 
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soft gluon and using Eqs. (Jg) and ([Tj)) we have 

( M^'\q, M) I M^'\q, M) ) + c.c. ^ - {gsfi'f 

M^'\{p}) I J(°)(g) ■ J'^W(g) I ^«(M) ) + c.c] (23) 

+ {gsfi'f [{M^'\{p})\J^^\q)-J^^'\q,e)\M^'\{p}))+ c.c. 

where c.c. denotes the complex conjugate. 

The first term on the right-hand side, which can be evaluated by using the expression 
(y) of the tree- level gluon current, has the same structure as Eq. (|T2|): 

m 

{M<^'\{p})\J^;\q)-J^'^'\q)\M^'\{p}))+ c.c. = 2 ^ S,,{q) |^g)(M)P , (24) 

where Sij{q) is the eikonal function in Eq. (|T3|). In particular, the two-particle colour 
correlations on the right-hand side are completely analogous to those at tree level (see 
Eqs. (|l^) and (|T^) ) and have been taken into account by defining the colour-correlated 
one-loop amplitude 



l-^S)(W)r = (-M^"HM)|T,-T,|A^W(|p}))+ e.c. . 



((0) 



(25) 



Using expression (|18|) for the one-loop contribution to the soft-gluon current, the second 
term on the right-hand side of Eq. (^) can be written as 



'M^'\{p}) I J^'\qyj^^'\q, e) | M^'\{p}) ) + c.c. 



1 {Any r3(l-e)r2(l + e) 



47r2 e2 



r(l-2e) 

• <! Ca cos(7re) [S^M]'^' l-MSS)(M)r (26) 



1,3 



+ 2sin(7re) Y^S,.{q) [S^MT (A., - A., - A,,) \Mfl,^{{p})\' 

where we have used {Xij — Xig — Xjqf = 1 for any possible configuration of incoming and 
outgoing momenta, and the notation stands for the sum over the different values of the 
indices {i ^ j,j ^ k,k ^ i). 



The one-loop contribution on the right-hand side of Eq. (^61) contains two terms. The 
first term only involves colour correlations between two hard partons, which are taken 
into account by the colour-correlated tree- amplitude A^|°j)({p}) defined in Eq. ([T^). Thus, 
apart from an overall e-dependent factor, its effect simply amounts to the rescaling Sij (g) — > 
[Sij{qy\'^'^ in the tree-level factorization formula (0). The second term instead has a 
different structure, because it leads to colour correlations between three different hard 
partons. These are included in the three-parton correlated tree-amphtude: 



i-^(L)(W)r = fa,c {M^'\{p}) I n I M^'\{p}) ) 



(27) 
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Note that the second term contributes only when there are four or more hard partons, 
because in the case of three partons colour conservation (see Eq. (^) can be used to show 
that the three-parton correlations vanish: 

fatcT^'T^T^ I M^'\puP2,P3) ) = -fatcT,^T^iT^ + T^) I A<(°)(Pi,P2,P3) ) = . (28) 

The absence of these correlations extremely simplifies the structure of the soft limit of the 
squared matrix elements with two or three hard partons. As a matter of fact, in these cases 
any product Tj ■ Tj can be expressed as a linear combination of Casimir operators (see 
Sect. H). Therefore, the colour algebra can exphcitly be carried out and the soft limit of 
\M.{q, is directly proportional to |A1({p})P at one-loop accuracy. The corresponding 

explicit expressions are presented in Sect. |^. 

Note also that in the limit e ^ the three-particle correlation contribution to Eq. (p^) 
is at most as singular as 1/e. More precisely, it gives rise to single poles 1/e only when there 
are two or more incoming partons. In fact, in the case of an outgoing soft gluon g, we have 
\j — \q — ^jq = +1 when i and j are both incoming momenta and Xij — Xiq — Xjq = — 1 
otherwise. Therefore, we can rewrite the second term in the curly bracket of Eq. (p6D as 

+4sin(7re) ' ^>^ril) i^MY fa^d M^'\{p}) l^T^T^ \ M^'\{p}) ) 

i,j{in),k 

-2sin(7re) Y^S^q) [SMY faUM^'\{p}) m^T; \ M^'\{p}) ) , (29) 

where the sum Yl'ij{m) k denotes the restriction of Yl!ijk ^'^ indices z, j of the incoming 
partons. Then, it is easy to show that the second term in Eq. (^9]) is of 0{e^) in the limit 
e ^ 0: 

-2sin(7re) J^S^q) [S^MY fabc{M^'\{p}) l^T^T^ \ M^'\{p}) ) = (30) 

-2vre 5^'5,,(g) faUM^'\{p}) \T^T,'T^ \ M^'\{p}) ) + 0{e') = 

+2ne J^SUq) faUM^"\{p}) \T^T^{T^^ + T^) \ M^'\{p}) ) + 0{e') = 0{e') .(31) 

i,k 

Here we have set [Sij{q)Y — > 1- Then we have used colour conservation, ^^^^^^T^ = 
-{T^ + T^), and the identity fabcT^Tf{T^ + T^) = for i ^ fc. 

4 Proof of factorization and 

calculation of the one-loop current 

In this section we derive the factorization formula (^) at one-loop order and we explicitly 
compute the one- loop contribution to the soft-gluon current. 

To simplify the analysis, it is convenient to work in a gauge with only physical gluon 
polarizations. We use the axial gauge n ■ A = with a light-like {v? = 0) gauge vector 
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n'^. The polarization vectors e^{k) of a gluon with momentum k thus fulfil the relations 
n-e{k) = and k-e{k) oc k'^. The sum over the gluon polarizations leads to the polarization 
tensor d^'^: 

d^^ik) = J2 ^'{k) e^{k) = ~g^^ + . (32) 

pol. 

The expression on the right-hand side corresponds to a dimensional-regularization scheme 
with — 2 = 2 — 2e gluon polarizations. For the sake of definiteness we use this explicit 
expression in all the intermediate steps of the calculation. However, we shall show that 
the final results are regularization-scheme independent by pointing out in which steps the 
scheme dependence might arise. 

4.1 Proof of factorization 

We study the soft behaviour of one-loop matrix elements by using the eikonal approximation 
and the soft-gluon insertion rules|] recalled above Eq. (p!0|). For this purpose, it is useful to 
decompose the one- loop matrix element MS^^ in three contributions, 

M^'^ = Mf^r,^ U^,,^ M%, , (33) 

which respectively represent the kinematic regions where the momentum k circulating in 
the loop is 

• hard: its momentum components are of the same order as those of the hard external 
momenta p, 

• collinear. k is parallel to one of the hard external momenta, 

• soft: its momentum components are much smaller than those of the hard external 
momenta p. 



For the purpose of the following discussion, we recall that the soft virtual behaviour of 
the one-loop matrix element A4^^\{p}) can be computed by using the soft-gluon insertion 
rules. The soft virtual gluon of momentum k is inserted (emitted and reabsorbed) on all 
the external legs (Fig. |l]) of the tree-level amplitude M.^'^\{p}), leading to the expression 



<yW)) = ^^?sV^^/|^^ [r,''\k)f r^^^\k) \M^'\{p})) , (34) 



where J^^^\k) is the (tree- level) soft-gluon current in Eq. (|]). 

"We use the eikonal approximation both for real and for virtual soft gluons. In the case of virtual 
gluons, the eikonal approximation is not always justified for each single Feynman diagram. Nonetheless, 
it is valid for any gauge-invariant set of Feynman diagrams, such as those computed in this section. This 
can be shown po| by simply using time-ordered perturbation theory, where the eikonal approximation is 
valid on a graph- by-graph basis. 
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Figure 1: Feynman diagrams that contribute to the soft behaviour of the one-loop amplitude. 
The shaded blob denotes the tree-level amplitude and the virtual gluon with soft momentum 
k can either (a) connect two different external legs i and j or (6) be emitted and reabsorbed 
by the same leg. 

We now consider the one-loop matrix element Ai^^\q, {p}) when the momentum q of 
the external gluon becomes soft. To apply the soft-gluon insertion rules, we perform the 
decomposition in Eq. (^) and discuss the three different kinematic regions of the virtual 
momentum k in turn. 

When k is in the hard region, all the internal lines in Ai^^\q, {p}) are highly off-shell. 
Thus the real soft gluon q can couple only to the external legs and we can factorize its 
contribution as in the case of tree-level amplitudes. Neglecting terms that are not singular 
in the soft limit g — 0, we get 

(«l M^Uq, W)) ^ 9sf^'e^{q)j;i^'\q) I ^LU(M)) , (35) 
where J^^^\q) is the soft-gluon current in Eq. (|]). 

We now consider the region where the loop momentum k is collinear to the momentum 
of one, say pi, of the hard external legs. Since we work in a physical gauge, the only 
diagram that is not dynamically suppressed [Q is that in which the loop leads to a self- 
energy contribution on the external leg pi (e.g. the diagram in Fig. |I] (6)). The soft gluon 
can then be inserted in this diagram in all possible ways: on the line pi before and after the 
self-energy contribution and on the self-energy lines themselves. However, we can exploit 
the colour coherence properties @] of QCD radiation. Since the loop-momentum k is parallel 
to Pi, the soft gluon q cannot distinguish the self-energy lines from a single line with total 
colour charge T". The sum of the insertions of q on this diagram is thus insensitive to the 
presence of the collinear loop momentum and leads (see, for instance. Sect. 3.4 in Ref. [0 
for a detailed similar discussion) to the same factor, T^Pi ■ s{q)/pi ■ q, as in the soft-gluon 
factorization at tree level. Considering the insertions of q on all the other external legs 
j 7^ i, we obtain a factorization formula of tree-level type also in the collinear region: 

(«l M^^J,{q, M)) ^ ^?5/iV'(g) J;(°)(?) I ^S.({P})) • (36) 

We finally have to deal with the region in which the loop momentum k is carried by a 
soft gluonQ. Unlike the case of the hard and collinear regions, where the one-loop effects 

** Quarks loops are dynamically suppressed when their momentum become soft. Without loss of infrared 
accuracy, we thus consider the quark loops as included in the hard or collinear regions. 
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can be factorized with respect to the tree- level current (g) (see Eqs. ( ^5|) and (^6l)), 
new 'non-factorizable' contributions appear when the loop momentum is soft. To single 
out these new contributions, we write the following identity: 

\M%M. M)) = 9s^^'e>^{q) Jf{q) \M%,{{p])) 

+ ( \M%M^ {P})) - 9s^^'e'^{<l)JfmM%m)) ) > (37) 

where we have added and subtracted the 'factorized' contribution. Then we combine the 
contributions from the hard, coUinear and soft regions by adding Eqs. (pSf), ( p6D and (pTf), 
and we obtain 

\M^'\q, M)) = 9s^^'e^{q) jf{q) \M^'\{p])) 

+ ( \M%M. {P})) - 9sf^'e^{q)jTi<l)\^Zi{p})) ) • (38) 

The first term on the right-hand side of Eq. (^) together with the contributions from 
Eqs. ([35| ) and (p6D have reconstructed the first term on the right-hand side of Eq. (|38D , 
which is exactly the first term on the right-hand side of the factorization formula (|17D . 
What remains to be done to prove the factorization formula is to relate the second term 
on the right-hand side of Eq. (0) with the contribution in the round bracket of Eq. (|38|) . 




Figure 2: Graphs that contribute to the one-loop soft current. 

For this purpose, we first note that when the real gluon q and the virtual gluon k are 
both soft, they can couple only to the external hard lines. In the corresponding Feynman 
diagrams, which are schematically represented by the first graph in Fig. H, the tree-level 
amplitude 7Vl''°-*({p}) is factorized in the soft limit. We can write: 

IM^tiqAp})) ^ i9sf^'f e^(g)K«(g,e) \M^'\{p})) , (39) 

where the kernel K^^^ (represented by the box in Fig. |^) denotes all the soft-gluon insertions 
of q and k on the hard-momentum lines. Then, we note that Ai^^\{p}) is factorized also 
in the expression (0) for ■M.^^Jft{{p}) ■ Therefore, the contribution in the round bracket 
of Eq. (^) can be recast in the form of the second term on the right-hand side of the 
factorization formula ([T^). Moreover, using Eqs. (|39j) and (^4]), we obtain the following 
explicit representation of the one- loop contribution J^^-* to the soft-gluon current (Fig. ^: 

.'^(g) J«(g,e)=.ng){K«(g,e)-j(f)(g)l [j'^\k)]' ■ J^^^\k)^ . 

(40) 
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4.2 Calculation of the one- loop current 



We now proceed to the explicit calculation of the one- loop soft current ' . Using the 
eikonal approximation, we have to evaluate the Feynman diagrams of the kernel K^^^ and 
to subtract those corresponding to the second term on the right-hand side of Eq. (^0|). We 
divide the diagrams in two classes: {A) the diagrams that involve interactions with a single 
hard line (Fig. and [B) all the remaining diagrams (Fig. ^). 




Figure 3: Feynman diagrams that depend on a single hard momentum. 

We first consider class {B). The diagrams (a), (6), (c) and [d) in Fig. | come from the 
kernel K'^^\ while Fig. ^ (e) represents the corresponding subtractions. We see that within 
this class there are diagrams that involve interactions between three different external lines 
i, J and /. These diagrams cancel. Indeed, the diagram in Fig. ^ {d) is exactly cancelled by 
that (e)-contribution in which q is emitted by the line / (the sum over the emissions of q is 
included in the factor J^^\q)). Thus the only non- vanishing terms coming from class (B) are 
those that involve interactions between two external lines, namely the diagrams (a), (6), (c) 
and the corresponding subtractions in (e). The diagrams (a), (c) and the subtractions in 
(e) are very similar: they can be combined in a simple way, because they only differ by the 
momentum and colour flow along the hard line i. The different factors coming from the 
line i in the diagrams (c) and (e) give 

' ' Piq + iO Pi{k + q) + iO ' ' piq + iO pik + iO ' ^ ' 

while the corresponding factor from the diagram (a) is 

_^J'brpa ^ I r^2) 

' ' Pik + iO pi{k + q)+iO ' 

Decomposing the colour factor in Eq. ( ^21) in its non-abelian and abelian components, 
T^T!^ = ifbacTi+T^Tj', and adding Eq. (p]), we see that the abehan component T^Tj' of the 
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Figure 4: Feynman diagrams that depend on two or three different hard momenta. 

diagram (a) exactly cancels the diagrams (c) and (e). In conclusion, the total contribution 
to the one-loop current from class (B) simply amounts to computing the diagram (6), 
which is non-abelian, and the non-abelian part of the diagram (a), which is obtained by 
the replacement TjTj'T"' —>■ ifbacTjT^ {h is the colour index of the virtual gluon) in its 
overall colour factor. 

The calculation of the diagrams in Fig. | (a) and (6) is quite simple (see below). Then 
we have to add the diagrams of class {A) (Fig. |^). An argument similar to that in Eqs. (p|) 
and ( ^2]) can be used to show that the abelian parts of these diagrams cancel^. Thus, only 
their non-abelian part has to be evaluated. Although straightforward, this calculation is 
quite cumbersome, in particular because it has to be carried out in the axial gauge n- A = 0. 
This cumbersome calculation can be short-circuited by exploiting gauge invariance. 

We first recall that the non-vanishing contributions from class {B) depend on the mo- 
menta and charges of (at most) two hard partons. Since the diagrams of class {A) involve 
interactions with a single hard line, we can split the one-loop soft current in two terms, 
JiP and J2P, that respectively denote the contributions that depend on one and two hard 
momenta: 

Jf^ (Q, e) = j;^lk<l, + j;^ik<l. • (43) 

All the diagrams of class (A) are included in Jip. Those of class (B) contribute to J2P 
and (because of possible cancellations in the dependence of one hard momentum) to Jip. 
We shall explicitly compute J2P and show that it is gauge-invariant. Thus Jip has to be 
gauge invariant as well, that is, it cannot depend on the gauge vector n^. Power counting 
and the fact that the one- loop current is non-abelian are then sufficient to determine Jip, 

^^More precisely, the abelian diagrams in Figs. || (a) and (6) are cancelled by the abelian part of the 
diagram in Fig. || (c). 
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apart from an overall factor /^^■'(e) that can only depends on e: 

m 

j;-?pi<l, e) = Ca f'\e) J2{p, ■ Tt . (44) 

i=i ' ^ 

This overall factor can then be obtained by imposing that the soft-gluon factorization 
formula be gauge-invariant (with respect to the polarizations of the real gluon q) or, equiv- 
alently, that the complete one-loop current in Eq. (^) be conserved. This leads to the 
constraint 

m 

q'jf\q. e) = q'j;Sk<l. ^) + Ca f'\e) Y^{p^ ■ g)"^ Tt = Q , (45) 

1=1 

which can be used to get f^^\e) (and, hence, J'^^lp) from the two-parton contribution J'^^2P- 
In particular, the calculation of the diagrams of class {A) can be avoided. 

We anticipate that our explicit calculation of the two-parton contribution J'^^2P gives 

g'" J^.2p(g, e) = 0. Thus the single-parton contribution J^^lp vanishes. Note that this does 
not mean that the diagrams of class (A) give a vanishing contribution, but rather that their 
contribution is cancelled by the terms of class (B) that depend on a single hard momentum. 

We can now complete our calculation of the one-loop current by explicitly computing 
the non-abelian part of the diagrams in Figs. ^ (a) and (6), and, more precisely, their 
contribution to </^.2^p. For the sake of definiteness, we write down the diagrams for the 
case in which all external momenta are outgoing. The diagram in Fig. ^ (a) depends on 
the polarization tensor d{k) of the virtual gluon. Using the expression in Eq. (^) and 
performing the Lorentz algebra of the spin numerator, the non-abelian part of the diagram 
in Fig. H (a) gives 



Pi ■ £{q) 



(A;2 + i(}){pj -k-iO) pt ■ {k + q) + iO 

Pi ■ Pj -Pj ■ k y Pi - k -J ) . (46) 



Pj-/c + iO\ k ■ n k ■ n 

The diagram^ in Fig. ^ (6) depends on the polarization tensors dik) and dik + q) of the 
virtual-gluon lines. Performing the Lorentz algebra of the spin numerator and exploiting 
the symmetry of the integrand under the transformation i ^ j, k ^ —{k + q), the diagram 
in Fig. I (6) gives 



■ P3 - P^P3 ■ (k + 2q) -prk (2P|:^ - Pf^^^^) + P^ik + qf 



{2nY (P + iQ)[p. . k - zO) pi-{k + q)+ iO {k + qf + zO 

Pi-n ^{k + 2q)-n\ ^ ^ Pj 



k ■ n 



(47) 



MNote that we use the eikonal approximation for the vertices and propagators of the hard-momentum 
hues i and j, while gluon propagators and the three-gluon vertex must be treated exactly. 
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By inspection of Eqs. (^) and (0), we see that the last term in the round bracket of Eq. (^) 
cancels the last term in the square bracket of Eq. (^T]). The remaining n-dependence 
of the integrands is proportional to pj ■ k. This factor cancels the fermion propagator 
l/{pj ■ k — iO), thus leaving a contribution that, although n-dependent, depends on a single 
hard momentum pi. This term can then be included, together with the class (A) diagrams, 
in the single-particle contribution j'^J^piq,^) to the one-loop current. We conclude that 

(the contribution of Eqs. (|^) and (|47|) to) the two-particle current jj^lp{q,e) is explicitly 
independent of the gauge vector n'^. Note also that ?7,-dependence cancels at the integrand 
level, and thus, the cancellation is completely insensitive to the actual prescription p2 



used to regularize the gauge pole l/{n- k) of the gluon polarization tensor (|3^) . 



Although the expression (|34 ), which we have used for the polarization tensor, corre- 
sponds to (i — 2 = 2 — 2e hehcity states for the gluon, our calculation also allows us to 
discuss the (in) dependence on the dimensional-regularization scheme. Other regularization 
schemes [pOj ^] use 2 helicity states for the gluon. At one-loop order, the difference even- 

in the integrands after having performed the Lorentz 

iBD and (^TD have no explicit 



tually amounts [^, ^ to set e 
algebra of the spin numerators. Since the integrands of Eqs. 



e-dependence, our result for e{q) ■ Jgp (and hence for the complete current e{q) ■ J''^^) does 
not depend on the scheme used to implement dimensional regularization. 

We proceed to the computation of e{q) ■ by adding the two-particle contributions 
from Eqs. (^Gf ) and (|47|). As for Eq. {^7\), the term proportional to k^ can be reduced to 
scalar integrals. Then using q ■ e{q) = and exploiting the symmetry with respect to the 
transformation i ^ j,k ^ 

^{q) ■ J2p\(l^^) = fabc 



— {k + q), we obtain 

d'^k Pi ■ e{q) 

Pi ■ Pj 



{2tiY pi-{k + q)+ iO (F + iO){pj -k-iO) 



+ 



1 



Pi ■ k + iO {k + q)'^ + iO \Pi ■ q 



Pi -Pj , o Pj'(l , 
k ■ q — Zpj ■ q Pi ■ k 



Pi ■ Q 



(48) 



The ultraviolet-finite integral can be easily performed in the coUinear frame in which pi and 
Pj are directed along the and '— ' light-cone directions, respectively. In this frame, the 
integration over the /c+-complex plane receives contributions from the pole l/{pj ■ k — iO) 
in the fermion propagators and from the poles in the gluon propagators. Closing the 
integration contour on the lower half-plane, and using the residue theorem, we can select 
only the gluon poles, which amounts to the following replacements in Eq. (0): 



1 



-2m 6+{k^ 



1 



Therefore Eq. (EHI) can be expressed as follows 



-2m6+{{k + qf 



(49) 



eiq)-r,j,'\q,e) = -'-U,yT,^T^P'-'^'^^ 



I{Pt^Pj\Pi\ 2p» ■ q) + IiPi,Pj;Pj; -2pj ■ q - iO) 



2 ^ ' Pi-q 
- I{q,Pj;Pj;-2pj ■ q - iO) - I{pj, q;pi;2pi ■ q) - I{pi, q;pj; -2pj ■ q - iO) 



(50) 
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in terms of the master integral 



I{p,p;r; s) 



X 0,2) P-P ^ 
{27iY-^ +^ ' {p-k){p-k) 2r-k + s 



(51) 



which depends on the hght-hke momenta p^,p^ and and on the complex scalar s. The 
evaluation of the d-dimensional integration in Eq. (|5lD over the on-shell gluon momentum 
k is straightforward and gives 



I(p,p]r] s) 



].n{p ■ r) {p ■ r) 
s"^ {p ■ p) 



-r2(l-e)r(l + e)r(l + 2e) . 



(52) 



Since I{p,p; r; s) vanishes when r = p or r = p, inserting Eq. ( [5^ ) in Eq. (|50|), we finally 
obtain 



IQn^ ^ ^ Pi-q \2pi-qpj- 



• ir2(i-e)r(i + e)r(i 



2e) 



IGvr 



1 + 



Pi 



'2m 



Pi 



Pi -Q Pj ■ <?/ V2pi ■ qpj -qe 



An Pi -pj 



1 T\l-e) r^(l + e) 
■ e2 r(l - 2e) 

where we have used the identity 

, , r(i + 6)r(i-e) 
= r(i + 2e)r(i-2e) 

and the symmetry under the replacement i ^ j. 

a(l), 



(53) 



(54) 



Note that the current in Eq. ( |55D is conserved, q^J^^2p{^i ^) — 0- Therefore, as discussed 
below Eq. (^5l), the single-particle term J^.^p vanishes and Eq. (|53|) gives the complete one- 
loop contribution to the soft-gluon current. We conclude that we have explicitly derived 
the result in Eq. ( [T8|) in the case of outgoing parton momenta PiiPj (when Ajj = A 



+1). The completely general result in Eq. (|T8|) is straightforwardly obtained by crossing 
{pi — >■ —pi and/or pj —Pj) the parton momenta in Eq. ( ^8|) and repeating the steps that 
lead to Eq. 



5 Processes with two and three partons 
and e+e"^ 3 jets at NNLO 

In general soft factorization formulae involve colour correlations (see Eqs. (|l^) and (p6D). 
As discussed at the end of Sect. |^, in the case of processes with two or three hard external 
partons, the correlations are completely given in terms of the products of the colour-charge 
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factors Tj ■ Tj. In this case, moreover, these colour-charge factors can be expressed in 
terms of the Casimir operators of the hard partons. Therefore the colour algebra can 
explicitly be carried out and we can obtain QED-like factorization formulae. We first 
recall the relations between colour-charge factors and Casimir operators and then give the 
factorization formulae. 

When there are only two hard partons in the amplitude, we can use colour conservation 
(Ti = -T2) to obtain 

Ti ■ T2 \Mipi,p2)) = -Tl \M{p,,p2)) = -Ci \M{p,,p2)) = -C2 \M{pi,p2)) , (55) 

where Ci = C2 is the Casimir of the hard partons. When the hard partons are three, the 
colour algebra can still be performed in closed form. Let us consider, for instance, the 
correlation term Ti ■ T2. Using colour conservation (T3 = — Ti — T2), its action on the 
ket \M{pi,p2,P3)) is 

2Ti ■ T2\M{pi,p2,P3)) = {Tl -Tl- Tl)\M{p,,P2,P3)) = (C3 - Ci - C2)\M{pi,p2,P3)) 

(56) 

and thus Ti ■ T2 is again given in terms of the Casimir invariants Ci. The same can be 
done for the other products Ti ■ T3 and T2 ■ T3. 

When the number of hard partons is four or more, colour correlations cannot be avoided. 
In fact the number of independent equations following from colour conservation is not 
sufficient to express the products Tj ■ Tj in terms of Casimir invariants. 

The colour-algebra results in Eqs. (|55| ) and (|56D can be inserted in Eqs. (0), (0) and 
(PB|) to straightforwardly obtain the soft limit of the squared amplitudes with a soft gluon 
and two or three hard partons (plus any number of colourless partons). In both cases, the 
contributions from Eqs. (|T^) and (^) reconstruct the one-loop expansion of the all-order 
squared amplitude |A^({p})P in Eq. (0), and we obtain the following factorization formulae 

\M{q,p,,p2)\' ^ glfi'^ACauiq) \M{p,,P2)\' , (57) 



2CiJi2(g) + C^(X23(g) +Xi3(g) - Ji2(g)) |A^(pi,P2,P3)I' 

(5^ 



where the kinematic factor Iij{q) includes the loop corrections to the eikonal function Sij{q) 
in Eq. ([T^) . To one- loop accuracy, we find 



J.,(g)=5.,(g)- 



gi 1 r^(l-e)r3(l+6) 



87r2 e2 r2(l -2e)r(l + 2e) 



(59) 



where we have used Eq. {p^. Note that in the amplitude with three hard partons, two 
of them (say, the ones with momenta pi and P2) have to form a particle-antiparticle pair 
(either a qq pair or a gluon pair) and the third one (with momentum p-^) has to be a gluon. 
In Eq. (|58|) we have therefore set Ci = C2 and C3 = Ca- 



The results in Eqs. (|57| ) and (0), when combined with the analogous QED-like factor- 
ization formulae for the emission of two soft gluons and of two 0, Q and three 
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1^ collinear partons, can be used to perform cross section calculations at NNLO for sev- 
eral processes. For instance, Eq. (0) is relevant to 2-jet production in e~^e~ annihilation 
and for the production of Drell-Yan, photon and vector-boson pairs in hadron collisions, 
while Eq. (^) can be used for the production of 3 jets in e~^e~ annihilation, (2+1) jets in 
deep-inelastic lepton-hadron collisions and (vector boson + jet) in hadron collisions. Ow- 
ing to the high precision of the data from e"'"e~-collider experiments, improved theoretical 
calculations for 3-jet production are highly demanded. The fact that colour and kinemat- 
ics factors are completely disentangled in the singular limits of the corresponding matrix 
elements certainly simplifies the structure of these calculations. 
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